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The spectral energy distribution (SED) of TeV blazars peaks both at keV 
and TeV energies. The X-ray emission is generally believed to originate in the 
synchrotron emission from relativistic electrons (and positrons) in the jet of these 
sources, while the origin of the gamma-ray emission is still being debated. We 

-vi | report results from a systematic study of X-ray spectral variability of Mrk 421 

and Mrk 501 during individual flares that last for several days, making use of 

£T} ', some of the high-quality data that have recently become available. The X-ray 

in 

their respective SEDs, so they together may offer additional insights into the 
O 



physical origin of X-ray variability. We modeled each of the time-resolved X- 
ray spectra over a single flare by adopting a homogeneous spatial distribution 
Q-f and an instantaneous power-law spectral distribution for the emitting particles. 

We focused on the variation of four key parameters: particle spectral index, 
maximum Lorentz factor, energy density, and magnetic field. While there is 
considerable degeneracy in the fits, we show that, in order to account for the 
X-ray spectral variability observed in Mrk 421, at least three of the parameters 
are required to vary in most cases, with the spectral index being one of them. 
The observations of Mrk 501 support the conclusion, although the quality of the 
data is not as good. We discuss the implications of the results. 

Subject headings: BL Lacertae objects: individual (Markarian 421 and Markar- 
ian 501) -- galaxies: active -- radiation mechanisms: non-thermal -- X-rays: 
galaxies 
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1. Introduction 

Blazars are known to be strong gamma-ray emitters. At TeV energies, BL Lac objects 
are the only type of blazars that have been detected so far. The spectral energy distribution 
(SED) of these TeV blazars invariably shows two characteristic peaks in the vF v represen- 
tation, with one located at X-ray energies and the other at TeV energies (e.g., Fossati et 
al. 1998). While the origin of the TeV emission is still not entirely clear, there is a general 
consensus that the X-ray emission is associated with the synchrotron radiation from highly 
relativistic electrons (and positrons) in the jet of blazars. Therefore, studying the X-ray 
properties of TeV blazars may, in a relatively model-independent manner, shed significant 
light on the properties of the emitting particles as well as of the environment that they are 
in. 

Blazars are highly variable at nearly all wavelengths. The variability is often in the form 
of discrete flares. Dramatic rapid flares have been observed of TeV blazars at X-ray energies 
(e.g., Cui 2004), as well as at TeV energies (e.g., Gaidos et al. 1996), with the former almost 
certainly being associated with the most energetic electrons in the jet. The observed X-ray 
flares span a wide range of timescales from many months down to less than an hour, and 
seem to form a well-organized hierarchy with weaker, shorter-duration flares superimposed 
on a stronger, longer-duration one (Cui 2004), perhaps indicating a scale-invariant physical 
origin (see also Xue & Cui 2005). The flaring phenomenon could be related to internal shocks 
in the jet (Rees 1978; Spada et al. 2001), or to the ejection of relativistic plasma into the 
jet (e.g., Bottcher et al. 1997; Mastichiadis & Kirk 1997), or to reconnection events in a 
magnetically dominated jet (Lyutikov 2003). 

While it is still difficult to pinpoint the exact mechanism at work, it may have become 
possible to investigate, quantitatively, which physical parameters evolve most strongly during 
a single flare on a timescale of several days, thanks to improvement in the quality of X-ray 
data in recent years. In this paper, we report results from such an investigation with two 
TeV blazars, Mrk 421 and Mrk 501, which are among the first sources discovered at TeV 
energies (Punch et al. 1992; Quinn et al. 1996; Bradbury et al. 1997). The sources were 
chosen for this work not only because of the availability of high-quality data but also of 
their diverse behaviors at X-ray energies. In particular, the observed X-ray spectra of the 
two sources fall on the opposite sides of the synchrotron peak in their respective SEDs. 
Modeling their spectra together may, therefore, help reveal common threads underlying the 
flaring phenomenon in TeV blazars. 
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2. Observational Aspects 

2.1. Data 

We observed Mrk 421 regularly with RXTE in 2003 and 2004, as part of a multiwave- 
length monitoring campaign (Blazejowski et al. 2005). The source became exceptionally 
bright in X-rays near the end of the campaign. Figure 1 shows the relevant portion of the 
(daily-averaged) ASM light curve. A number of distinct flares are present, with the strongest 
one peaking roughly at about 10 c/s (or 135 mcrab in the ASM band of 1.5-12 keV)! For 
detailed spectral studies, we used data from the pointed RXTE observations taken with the 
Proportional Counter Array (PCA). The times of these observations are noted in Fig. 1. The 
PCA detector consists of five nearly identical proportional counter units (PCUs). However, 
only PCU and PCU 2 were in use throughout the campaign. PCU has lost its front 
veto layer, so the data from it are prone to contamination by events caused by low-energy 
electrons entering the detector. For this work, therefore, we chose PCU 2 as our "standard" 
detector for flux normalization and spectral analysis. We also only used data taken in the 
Standar d2-mode, which have a time resolution of 16 seconds and cover the full energy range 
(nominally 2-60 keV). 

We followed Cui (2004) closely in reducing and analyzing the PCA data. Briefly, the data 
were reduced with ftools v5.2. For a given observation, we first filtered data by following the 
standard procedure for faint sources, 1 which resulted in a list of good time intervals (GTIs). 
We then simulated background events for the observation by using the latest background 
model (pca_bkgd_cmfaintl7_eMv20031123.mdl) that is appropriate for faint sources. Using 
the GTIs, we proceeded to make a spectrum for each PCU by using data from only the 
first xenon layer (which is most accurately calibrated), which limits the spectral coverage 
to roughly 3-25 keV. Since few counts are detected at higher energies, the impact of the 
reduced spectral coverage is very minimal. We repeated the steps to derive a corresponding 
background spectrum for the PCU from the simulated events. 

Unlike Mrk 421, Mrk 501 has been fairly quiet in recent years. For this work, we obtained 
archival RXTE data taken during the 1997 outburst of the source. For direct comparisons 
we treated the data in the exact same manner as in the case of Mrk 421. 



^ec the online RXTE Cook Book at http://heasarc.gsfc.nasa.gov/docs/xte/recipes/cook_book.html. 
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2.2. Results 

Figure 2 shows the PCA light curves of Mrk 421 and Mrk 501, covering the respective 
periods of strong flares. The coverage gaps in both cases are mainly caused by the require- 
ment for coordination with ground-based facilities. For this work, we would ideally like to 
focus on one single flare and follow the spectral evolution of the source through both the 
rising and decaying periods. We were able to do so in the case of Mrk 421. A large flare was 
well sampled by the observations, as shown in Fig. 2, and we chose five of the observations 
for further analyses. As pointed out by Cui (2004), however, it is, strictly speaking, not 
possible to isolate a flare from the flaring hierarchy. The notion of "single" flares is therefore 
meaningful only when a timescale is specified. Here, we are specifically interested in flares 
that last for days. In the case of Mrk 501, no single flare was adequately sampled. We se- 
lected representative observations that covered the rising phase of one flare and the decaying 
phase of a subsequent flare, as also shown in Fig. 2. Note that Mrk 501 was much fainter 
than Mrk 421. 

For each observation, we jointly fitted the PCU spectra in XSPEC vll.3.1 (Arnaud 
1996). An additional multiplicative factor was introduced to take into account any uncali- 
brated difference in the overall throughput of the PCUs, with reference to PCU 2. Following 
the usual practice, we added 1% systematic uncertainty to the data. Given the lack of sen- 
sitivity of the PCA below about 2 keV, we adopted the Galactic values for the line-of-sight 
hydrogen column densities (1.38 and 1.8 x 10 20 cm" 2 for Mrk 421 and Mrk 501, respectively; 
Dickey & Lockman 1990) and fixed them in the fitting process. 

We experimented with three empirical models: power law, broken power law, and power 
law with an exponential cut-off. The cut-off power law provides best fits to the data for the 
selected observations of Mrk 421 (as shown in Fig. 3), although fits with the broken power 
law are also acceptable. For Mrk 501, the broken and cut-off power law models fit the data 
almost equally well. In both cases, we used the best-fit model to unfold each observed (count) 
spectrum to derive the corresponding SED for further theoretical modeling. We should note 
that only the PCU2 data were used to produce the SED, while data from all active PCUs 
were used to find the best-fit model. 

Figure 3 shows representative X-ray SEDs for Mrk 421 and Mrk 501, respectively. Both 
sources exhibit significant spectral variability during the flares. In general, the spectrum 
seems to harden toward high fluxes. It is also clear, from the figure, that the measured X- 
ray SEDs of Mrk 421 and Mrk 501 lie on the opposite sides of their respective (synchrotron) 
peaks, which are (mostly) below 2 keV in the former case but above 20 keV in the latter. 
Therefore, together, the sources should be representative of TeV blazars. 
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3. Theoretical Aspects 



To gain insights into variability of the underlying distribution of electrons, as well as 
their environment, we proceeded to construct a physical model and to apply the model to 
the data in a systematic manner. 



3.1. Model 

We assume that each flare of interest originates in one localized region in the jet and 
that the electron and magnetic field are distributed homogeneously in the region. We further 
assume that the energy spectrum of emitting electrons is of power-law shape, 

71(7)^7 = iVo7~ p <i7, 7 min < 7 < 7 max . (1) 

The total energy density of electrons is then 

E tot = AT m e c 2 7minP ~^ , p^2. (2) 

p — 2 

For an electron with Lorentz factor of 7, the differential power of its synchrotron radi- 
ation is given by 



dP{u) 27T V / 3e 2 z/ L 



dv 



- / K 5/3 (t)dt 

V c J v/v c 



ergs-'Hz" 1 ), (3) 



where the Larmor frequency v L = j--^-, the critical frequency v c = |^l7 2 , and K n (t) is the 
modified Bessel function of the second kind of order n. Assuming that the flaring region is 
of spherical shape with radius r, the spectrum of the integrated emission (in the jet frame) 
is 

4 <■ T max , sdPUA , 

^ 7min 

If the Doppler factor of the jet is S, the observed spectrum is given by 

FobsK bs ) = 5 3 F„(v ohs /5) (5) 

It is often taken casually that such an optically-thin synchrotron spectrum is also of power- 
law shape with a spectral index of a — (p — l)/2. This is only true when the portion of the 
spectrum of interest is far away from the critical frequencies corresponding to either 7 m i n or 
7max- The issue is of particular relevance to the case of Mrk 421 because its X-ray spectrum 
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seems to be entirely associated with electrons near 7 max , as shown in Fig. 3. Therefore, care 
must be taken in the calculation. 

There are many parameters in the model. The objective of this work is to systematically 
investigate the roles of "key" parameters in producing the observed X-ray spectral variability 
during a major flare. We identified these parameters as: p, 7 max , E tot , and B. We adopted 
E tot rather than Nq because (1) it seems obvious that iVo must vary during a flare, so few 
insights would be gained; and (2) E tot , though depending on other parameters of interest, 
is a well-defined physical quantity and provide important information on the energetics of 
the flaring process. All other parameters were fixed during the fitting process. The size of 
the emitting region was chosen to be compatible with the observed timescale of the flare, 
r = 1 x 10 16 and 5 x 10 15 cm for Mrk 421 and Mrk 501, respectively; the Doppler factor 
was set at a nominal value for TeV blazars (5 = 15); and the minimum Lorentz factor of the 
electrons was made sufficiently low (7 m i n = 10 4 ) such that it does not affect the conclusions 
reached. For this work, we adopted the luminosity distances of d = 137 and 156 Mpc for 
Mrk 421 (z = 0.030) and Mrk 501 [z = 0.034), respectively, which were derived with the 
following cosmological parameters: H = 65 km s _1 Mpc -1 , Q m = 0.3, and fix = 0.7. 



3.2. Results 

For each spectrum shown in Fig. 3, a grid search was carried out to find statistically 
acceptable solutions. A solution is defined as a combination of p, 7 max , -Etot, and B that 



leads to a fit with reduced x 2 m the range of 1 ± \j2/v, where v is the degree of freedom of 
the fit. Note that a x 2 distribution can be approximated by a normal distribution with mean 
of v and standard deviation of sfh) when v is sufficiently large. Here, v is typically 30-40. 
The ranges of the parameters covered in the search are: p in 1.00-4.00, B in 10~ 3 -10 2 G, 
7 max in 10 5 -10 11 , and E tot in 10 8 -10 15 cm -3 . The step used is linear for p, with Ap = 0.01, 
but logarithmic for B, 7 max , and E tot , with 50 and 25 points per decade for Mrk 421 and 
Mrk 501, respectively. 

Table 1 summarizes the results for Mrk 421 and Mrk 501 separately. The ranges shown 
indicate the portion of the parameter space where solutions are found. It is quite clear that 
p is quite well constrained in both cases. For illustration, Fig. 4 shows the distributions of 
p for Mrk 421. The p distribution is nicely peaked. The most probable values of p are also 
shown in Table 1 for both sources. The results seems to indicate that the electron spectrum 
tends to harden toward the peak of a flare. There is a substantial amount of degeneracy 
among the other three parameters: 7 max is fairly well constrained in Mrk 421 but not at all 
in Mrk 501 (presumably due to poorer data quality); B is poorly constrained in both cases; 



-7- 



and E tot is not constrained at all. The distributions of these parameters are fairly flat with 
no dominant peaks. The best fits are shown in Fig. 3. 

To break the degeneracy, we proceeded to examine specialized cases, in which a subset 
of the parameters are fixed at some "typical" values, as seen from modeling broadband SEDs 
of the sources. It became immediately clear that, to account for the observed X-ray spectral 
variability of Mrk 421, p cannot be fixed. This is not the case for Mrk 501, presumably due 
to the poorer statistics again. Taking the view that the same processes are operating in both 
sources, we allowed p to vary in all subsequent analyses. 

We first experimented with fixing one of the remaining parameters, B, 7 max , and Etot- 
We found acceptable fits to all of the spectra in all cases. As examples, Table 2 and Table 3 
show the results obtained with B = 0.10, 0.83 G and 7 max = 5.0 x 10 5 , 5.2 x 10 6 for Mrk 421 
and Mrk 501, respectively (see, e.g., Blazejowski et al 2005 and Pian et al. 1998). The free 
parameters are now much more constrained. We then proceeded to investigate whether we 
could fix two of the parameters and still obtain solutions. We found that it was not possible 
to account for the observed spectral variability in Mrk 421, except for the case in which B 
and 7 ma x were fixed, which confirms that the synchrotron peak of the SED (oc 7m ax -^) did 
not vary too much during the flare. As for Mrk 501, the data are much less constraining; we 
were still able to find satisfactory fits to the spectra. For completeness, we also examined 
the possibility of fixing all three parameters. In this case, we failed to obtain acceptable fits 
to all of the spectra for either source. 

It is often assumed that the ratio of the energy of the radiating particles to that of the 
magnetic field remains constant, in theoretical studies of gamma ray bursts and blazars. To 
check if the assumption can be accomodated by the data, we looked at cases in which E to t/B 2 
remains constant. For both sources, acceptable solutions were found. As an example, we 
show in Table 4 those solutions selected from the full set (as shown in Table 1) that have 
Etot/B 2 fixed at some fiducial values. In this case, the other parameters are even more 
constrained, although it is not possible to know the exact value of the ratio a priori. 



4. Discussion 

We have obtained time-resolved X-ray spectra of Mrk 421 during the rise and fall of a 
prominent X-ray flare, as well as of Mrk 501 during two consecutive flares (on a timescale of 
a few days) for comparison. In general, the X-ray spectrum steepens as the flux decreases 
(see Fig. 3), as was seen previously in X-ray selected BL Lac objects in general (Giommi 
et al. 1990). For Mrk 421, the synchrotron peak is clearly below the low-energy sensitivity 



threshold of the PCA detector (~2 keV), judging from the shape of the SEDs, except for 
one case ("£3" in Fig. 3) in which the peak appears to have moved into the PCA passing 
band. A shift in the synchrotron peak with flux has been seen in Mrk 501 previously (Pian 
et al. 1998; Sambruna et al. 2000). In our case, the synchrotron peak of Mrk 501 is clearly 
above the spectral range of the data shown in Fig. 3 (but not necessarily above that of the 
PCA), which is known (e.g., Pian et al. 1998; Catanese et al. 1997; Lamer & Wagner 1998), 
and is also why we selected it for comparison with Mrk 421 in the first place. Because of 
that, our SED modeling suffers from serious degeneracy between B and 7 max . For instance, 
the unusually- high values allowed for 7 max or unusually-low values for B is the consequence 
of such uncertainty. 

To cast light on the origin of the observed X-ray spectral variability of both Mrk 421 and 
Mrk 501, we took an approach to model instantaneous properties of the sources throughout 
the flares. In other words, we were not concerned with the steady-state properties. For 
simplicity, we assumed that the emitting electrons always follow a power-law spectral dis- 
tribution, which may or may not be the case in reality, after heating and cooling processes 
are properly taken into account. Moreover, we assumed that the spatial distribution of the 
electrons is homogeneous during a flare. Even with such an idealistic approach, the model 
still contains too many parameters to de-couple completely. We narrowed our focus onto 
the evolution of four key parameters: p, 7 max , E tot , and B, with other parameters fixed at 
nominal values derived from variability timescales and broadband SEDs (in particular, the 
TeV gamma-ray properties). In spite of serious degeneracy in the solutions, we are able to 
reach a few conclusions that we think are robust. 

As expected, the high-quality data of Mrk 421 have provided the most severe constraints 
on the model. To account for the observed spectral variability of the source, we found that 
the electron spectral index (p) must vary. There appears to be a general trend of decreasing p 
toward high fluxes. The only exception seen (t x in Fig. 4) is probably due to "contamination" 
by a previous flare (see Fig. 2). Of course, it could also be associated with local variability 
within the flare of interest. Again, there may always be smaller flares superimposed on it 
(Cui 2004). The change in p alone is, however, not sufficient for explaining the observed 
X-ray spectral variability. At least two of the three remaining parameters must also vary, in 
all but one combinations. The lone exception is that acceptable solutions do exist even when 
B and 7 max are both fixed. This is an indication that the synchrotron peak did not vary too 
much during the flare. Even when all four parameters are allowed to vary, both p and 7 max 
are reasonably well constrained (see Table 1). If we set B or 7 max at typical values from 
modeling the broadband SED, as well as time variability, in the context of the synchrotron 
self-Compton scenario, the other parameters are much more constrained and their values are 
also quite typical (see Tables 2 and 3). 
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The quality of data is not nearly as good for Mrk 501 as for Mrk 421. However, we 
feel that it is a worthwhile exercise to see whether the same model can also be applied 
to a different TeV blazar, like Mrk 501 whose observed X-ray spectrum is almost entirely 
below the synchrotron peak. We followed precisely the same procedures in modeling the 
X-ray spectra of Mrk 501 as Mrk 421. From Tables 1, 2, and 3, it is clear that the model 
parameters are poorly constrained in this case. One interesting (and solid) result is that 
the electron spectrum is much harder for Mrk 501 than for Mrk 421, even after taking into 
account the large uncertainties. Moreover, the electron spectrum hardens toward high fluxes 
during the rising phase of the first flare, although the trend is not as clear during the decaying 
phase of the second flare, perhaps due to the presence of more complicated local variability 
(see Fig. 2). Because the synchrotron peak is beyond the energy range covered by the data 
in this case, B and 7 max are hardly constrained. 

We should stress that it is not our goal to exhaustively explore the entire parameter space 
of the model. The scope of our investigation is limited only to a few selected parameters, 
which could affect the conclusions. First of all, although the effects of some of the parameters 
not selected (e.g., the size of the emitting region) can be absorbed into E to t, those of others 
may still contribute to the observed X-ray spectral variability. For instance, the Doppler 
factor of the jet might not remain constant during a flare and thus lead to not only variable 
Doppler boosting of the source intensity (which only affects the overall normalization) but 
also variable Doppler shift of the photon energy (see Eq. 5). We should note, however, 
that we also ran a case for Mrk 421 where the size of the emitting region and the Doppler 
factor were fixed at very different values (R = 10 15 cm and S = 50). Although the exact 
values of the best-fit parameters are different, our conclusions remain the same. Second, 
the homogeneous electron distribution considered here may be a gross over-simplification 
of the reality. Multiple electron populations may be present, with different properties and 
environments, and all may contribute to the X-ray emission. However, in our view, it seems 
more likely that a single rapid flare is associated with one electron population, which is why 
we tried to focus on specific flares. The most important conclusion of this study is that 
the observed X-ray spectral variability of TeV blazars is likely caused by a combination of 
spectral variation of the radiating electrons and changes in the environment (e.g., B). This 
is hardly surprising because it seems unlikely that, in reality, the triggering of a flare involves 
changes only in one of the parameters. Moreover, our results show that the spectral index (p) 
of the particle distribution must vary during a flare, together with the overall normalization 
(No). 

This work was supported in part by NASA Grants NAG5- 13736 and NNG04GP86G. 
It has made use of data obtained through the High Energy Astrophysics Science Archive 
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Table 1: Spectral Solutions: Full Set 








Observation MJD p 
(Data Set) 


B 

(Gauss) 


Tmax 


Etot 
(cm -3 )" 


mpv(min,max) 6 


(min,max) 


(min,max) 


(min,max) 



421 £1 

421 n 

421 £3 
421 £4 
421 £5 



53114.25 
53115.31 
53116.30 
53118.24 
53119.25 



3.16(3.04,3.27) 
3.45(3.20,3.64) 
2.76(2.54,2.93) 
3.39(3.13,3.57) 
3.87(3.26,3.95) 



0.002,3.98 
0.003,5.25 
0.002,4.57 
0.003,5.50 
0.002,4.37 



1.9xl0 5 ,8.3xl0 6 ) 
1.5xl0 5 ,6.0xl0 6 ) 
1.7xl0 5 ,9.5xl0 6 ) 
1.4xl0 5 ,6.0xl0 6 ) 
1.7xl0 5 ,1.0xlO n ) 



;io 8 ,io 15 ) 
;io 8 ,io 15 ) 
;io 8 ,io 15 ) 
;io 8 ,io 15 ) 

10 8 ,10 15 ) 



501 Tl 
501 T2 
501 T3 
501 T4 
501 T5 
501 T6 
501 T7 



50553.37 

50554.17 
50554.44 
50580.36 
50581.19 
50581.36 
50582.36 



2.60(1.77,2.73) 
2.28(1.35,2.46) 
2.13(1.04,2.45) 
2.56(1.76,2.71) 
2.30(1.39,2.71) 
2.07(1.57,2.37) 
2.63(1.32,2.72) 



0.002,75.9 
0.002,100. 
0.002,100. 
0.002,91.2 
0.002,83.2 
0.002,6.31 
0.002,75.9 



1.7xl0 5 ,1.0xlO n ) 
1.7xl0 5 ,1.0xlO n ) 
1.6xl0 5 ,1.0xlO n ) 
1.7xl0 5 ,1.0xlO n ) 
1.6xl0 5 ,1.0xlO n ) 
1.9xl0 5 ,1.3xl0 7 ) 
1.6xl0 5 ,1.0xlO n ) 



(10 8 ,10 15 ) 
(10 8 ,10 15 ) 
(10 8 ,10 15 ) 
(10 8 ,10 15 ) 
(10 8 ,10 15 ) 
(10 8 ,10 15 ) 
(10 8 ,10 15 ) 



a E tot has been scaled by the electron rest energy 
b mpv: most probable value 



Table 2: Spectral Solutions: Subset with Fixed B 



Observation 
(Data Set) 



MJD 



P 



B 

(Gauss) 



7i. 



E< 



tot 



cm 







mpv(min,max) 


(fixed) 


(min,max) 


(min,max) 


421 £1 


53114.25 


3.11(3.05,3.23) 


0.10 


(I.lxl0 6 ,1.3xl0 6 ) 


(2.5xlO n ,3.5xlO n ) 


421 £2 


53115.31 


3.48(3.22,3.62) 


0.10 


(1.0xl0 6 ,1.3xl0 6 ) 


(4.6xlO n ,l.lxl0 12 ) 


421 £3 


53116.30 


2.75(2.62,2.89) 


0.10 


(I.lxl0 6 ,1.3xl0 6 ) 


(2.3xlO n ,3.2xlO n ) 


421 £4 


53118.24 


3.34(3.14,3.56) 


0.10 


(1.0xl0 6 ,1.3xl0 6 ) 


(4.6xlO n ,l.lxl0 12 ) 


421 £5 


53119.25 


3.88(3.30,3.95) 


0.10 


(9.5xl0 5 ,1.0xlO n ) 


(2.5xlO n ,1.2xl0 12 ) 


501 Tl 


50553.37 


2.60(1.91,2.71) 


0.83 


(4.8xl0 5 ,1.0xlO n ) 


(4.8xl0 9 ,2.3xl0 12 ) 


501 T2 


50554.17 


2.28(1.75,2.34) 


0.83 


(6.3xl0 5 ,1.0xlO n ) 


(8.3xl0 9 ,7.6xl0 13 ) 


501 T3 


50554.44 


1.97(1.04,2.31) 


0.83 


(4.4xl0 5 ,2.5xl0 6 ) 


(6.9xl0 9 ,2.3xl0 10 ) 


501 T4 


50580.36 


2.56(1.94,2.63) 


0.83 


(5.2xl0 5 ,1.0xlO n ) 


(6.3xl0 9 ,2.8xl0 12 ) 


501 T5 


50581.19 


2.42(1.62,2.63) 


0.83 


(4.4xl0 5 ,5.8xl0 9 ) 


(5.2xl0 9 ,2.8xlO n ) 


501 T6 


50581.36 


2.04(1.78,2.31) 


0.83 


(4.8xl0 5 ,6.9xl0 5 ) 


(5.2xl0 9 ,6.9xl0 9 ) 


501 T7 


50582.36 


2.63(1.38,2.67) 


0.83 


(3.6xl0 5 ,9.1xl0 10 ) 


(3.6xl0 9 ,6.3xlO n ) 
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Table 3: Spectral Solutions: Subset with Fixed 7 n 



Observation 
(Data Set) 



MJD 



P 



B 

(Gauss) 



7u 



Etot 
fcm _3> ! 



mpv(min,max) 


^min,max) 


(fixed) 


(min,max) 


421 £1 


53114.25 


3.09(3.07,3.12) 


;o.52,o.58; 


5.0xl0 5 


(6.6xl0 9 ,7.6xl0 9 ) 


421 £2 


53115.31 


3.45(3.23,3.60) 


;o.4o,o.69; 


5.0xl0 5 


(I.lxl0 10 ,2.5xl0 10 ) 


421 £3 


53116.30 


2.76(2.59,2.89) 


;o.5o,o.76; 


5.0xl0 5 


(4.6xl0 9 ,8.3xl0 9 ) 


421 £4 


53118.24 


3.39(3.13,3.56) 


;o.4o,o.69; 


5.0xl0 5 


(1.3xl0 10 ,2.4xl0 10 ) 


421 £5 


53119.25 


3.80(3.29,3.92) 


;o.36,4.i7; 


5.0xl0 5 


(1.0xl0 8 ,1.6xl0 10 ) 


501 Tl 


50553.37 


2.58(1.88,2.70) 


;o.oi,i5.8; 


5.2xl0 6 


(1.0xl0 8 ,8.3xl0 13 ) 


501 T2 


50554.17 


2.26(1.35,2.34) 


;o.oi,3o.2; 


5.2xl0 6 


(1.0xl0 8 ,1.0xl0 14 ) 


501 T3 


50554.44 


1.98(1.04,2.32) 


;o.oi,o.25; 


5.2xl0 6 


(2.8xlO n ,1.4xl0 14 ) 


501 T4 


50580.36 


2.55(1.76,2.62) 


;o.oi,i7.4; 


5.2xl0 6 


(1.0xl0 8 ,1.0xl0 14 ) 


501 T5 


50581.19 


2.29(1.60,2.58) 


;o.oi,o.o5; 


5.2xl0 6 


(2.8xl0 12 ,1.3xl0 14 ) 


501 T6 


50581.36 


2.04(1.58,2.36) 


;o.oi,o.o2; 


5.2xl0 6 


(2.1xl0 13 ,1.2xl0 14 ) 


501 T7 


50582.36 


2.60(1.32,2.66) 


;o.oi,i3.2; 


5.2xl0 6 


(1.0xl0 8 ,1.6xl0 14 ) 



Table 4: Spectral Solutions: Subset with Fixed E tot /B 2 



Observation 
(Data Set) 



MJD 



P 



B 

(Gauss) 



7u 



K 



tot 



cm 



mpv(min,maxj 



mm, max) 



Jmm,max) 



mm 



max 



421 £1 

421 £2 
421 £3 
421 £4 
421 £5 



53114.25 
53115.31 
53116.30 
53118.24 
53119.25 



3.11(3.05,3.16) 
3.49(3.31,3.59) 
2.78(2.70,2.86) 
3.36(3.23,3.52) 
3.83(3.30,3.91) 



0.10,0.11) 
0.12,0.14) 
0.10,0.11) 
0.12,0.14) 
0.10,0.14) 



I.lxl0 6 ,1.2xl0 6 ) 
9.5xl0 5 ,l.lxl0 6 ) 
1.2xl0 6 ,1.3xl0 6 ) 
9.5xl0 5 ,1.0xl0 6 ) 
9.5xl0 5 ,1.0xlO m 



I.lxl0 6 ,1.0xl0 n ) 
1.9xl0 6 ,1.0xlO n ) 
;i.0xl0 6 ,2.8xl0 6 ) 
9.1xl0 5 ,1.0xlO n ) 
9.1xl0 5 ,5.2xl0 8 ) 
1.0xl0 6 ,1.0xl0 6 ) 
S.3xl0 5 ,3.3xl0 10 ) 



2.5xlO n 
3.6xlO n 
2.5xlO n 
3.6xlO n 
2.5xlO n 



2.8xlO n ) 
4.8xlO m 
2.8xlO m 
4.8xlO n ' 
4.8xlO n ' 



501 Tl 
501 T2 
501 T3 
501 T4 
501 T5 
501 T6 
501 T7 



50553.37 
50554.17 
50554.44 
50580.36 
50581.19 
50581.36 
50582.36 



2.60(2.27,2.71) 
2.28(2.15,2.34) 
1.85(1.70,2.27) 
2.56(1.92,2.62) 
2.46(2.05,2.60) 
2.18(2.15,2.22) 
2.63(1.82,2.66) 



0.28,1.20) 
0.33,3.31) 
0.30,0.36) 
0.28,1.20) 
0.28,0.58) 
0.28,0.28) 
0.25,0.76) 



5.8xl0 10 
8.3xl0 10 
6.9xl0 10 
5.8xl0 10 
5.8xl0 10 
5.8xl0 10 
4.8xl0 10 



l.lxlO 12 ' 
8.3xl0 12 ' 
l.OxlO 11 ' 

l.lxlO 12 ' 
2.5xlO n ' 
5.8xl0 10 ) 
4.4xlO n ' 
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Fig. 1. — Daily-binned ASM light curve of Mrk 421 in 2004. Each data point represents 
a weighted average of raw count rates in the summed band (1.5-12 keV) over 1 day. For 
clarity, the error bars are not shown. The times of the pointed PCA/ RXTE observations 
are indicated with vertical lines at the top. The arrow points to a flare that was very well 
covered by the PCA observations. For reference, the Crab Nebula produces about 75 c/s. 
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Fig. 2.— PCA light curves of Mrk 421 in 2004 (left) and Mrk 501 in 1997 (right). The 
results were derived from the PCU 2 data alone. The count rates were computed in the 2-60 
keV band. For reference, the Crab Nebula produces about 2700 c/s/PCU. The observations 
selected for subsequent analyses are marked. Note the presence of coverage gaps. 
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Fig. 3. — Representative X-ray spectra of Mrk 421 (left) and Mrk 501 (right). The solid 
lines show the best-fit models. The parameters are: for Mrk 421, (tl: p = 3.16, B = 0.22, 
7max = 8.2 x 10 5 , E tot = 1.05 x 10 11 ), (t2: p = 3.45, B = 0.25, 7max = 7.1 x 10 5 , E tot = 
1.75 x 10 11 ), (tS: p = 2.72, B = 0.21, 7max = 8.3 x 10 5 , E tot = 1.05 x 10 11 ), (t£ p = 3.38, 
B = 0.31, 7 max = 6.3 x 10 5 , E tot = 1.1 x 10 11 ), (*5: p = 3.71, 5 = 0.29, 7 max = 7.9 x 10 5 , 
E tot = 0.95 xlO 11 ); for Mrk 501, (Tl: p = 2.52, B = 0.10, 7max = 3.1 x 10 6 , E tot = 5.0xlO n ), 
(T3: p = 2.07, B = 0.10, 7max = 2.2 x 10 6 , E tot = 6.5 x 10 11 ), (TJ,: p = 2.44, B = 0.11, 
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Fig. 4. — Normalized distributions of the electron spectral index for Mrk 421. The results 
shown were derived for the case in which all key parameters are allowed to vary (see text). 



